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Abstract 



> 

o 

\& \ By means of the inverse techniques we analyse the evolution of purely spin-1/2 systems 

in homogeneous magnètic fields as well as the gener ation of exact solutions. Some evolution 
loops, dynamical processes for which any state evolves cyclically, are presented, and their 
í> | corresponding geomètric phases are evaluated. 

q ■ In the direct approach to quantum dynamics, the main goal is to look for solutions \ifi(t)) G 

\0 ', S(TC) to Schròdinger's equation 

g;: ihj t m)) = H(t)\m), (i) 

provided the Hamiltonian H(t) and the initial state |t/>(0)) are given, where S(7i) = {\ip) G 
TC | (i/)\ip) = 1}, Tï is the system Hilbert space [1]. The same ideology is present in the study of 



the 'shadow' motion of \ip(t)) on projective Hilbert space V, where p G V represents a possible 
physical state. V is obtained through a projective map ix : S(TÍ) — > V, and all = e^\ip), (fi G R 
q-i (a ray of S(TC)), are carried into the same p = 7r(\if))) G V. In this way, S(TC) is seen as a fibre 
bundle with base space V and fibre U{\) [2-3]. 

For spin-1/2 systems in magnètic fields B(t), the spin part of the Hamiltonian reads: 

H(t) = -l·iB(t) ■ S = -^-B(í) • cr, (2) 

where S is the spin operator, cr-components are Pauli matrices, S(7í) is a sphere S 3 embeded in 
R 4 , V is a sphere S 2 embeded in R 3 , and ^i\ip)) = n = (ifi\cr\i/;) G S 2 . The evolution equation 
for n(t) is obtained from (1-2): 

^ = -pB(t) x n(í). (3) 

Equations (1-3) are usually solved by means of numerical integration. 

It is interesting also to look for exactly solvable models where the relevant information would 
be found in one analytic expression. However, it is hard to find them due to the complexity to 
sum up the 'continuous' Baker-Campbell-Haussdorf exponent arising when the evolution operator 
is a true exponential [4]. It is well known that exact solutions to (1-3) arise when: i) the field has 
a fixed direction, B(í) = B(t)eB] ü) B(í) rotates around a fixed direction with constant angular 



velocity, and |B(í)| = |B(0)|. In the first case, by solving the eigenproblem for the constant 
matrix ob = cr ■ e# and expanding then \tp(0)) in that basis the general solution \ip(t)) is found. 
For the rotating field, the key is the 'transition to the rotating frame' where it arises a constant 
effective Hamiltonian, which can be solved as previously [5]; by coming back to the initial frame 
the solution is gotten. 

For general B(í) one can try again the transition to the rotating frame in which the direction 
of B(í) would be fixed though the effective Hamiltonian won't be constant; thus, this problem 
generally cannot be solved (see, e.g., [6]). This is the motivation to look for alternative techniques 
to circumvent these difficulties. 

In this work, a method to generate solutions for (1-3) is studied. We will use inverse techniques, 
i.e., we are going to look for the magnètic fields inducing a given n(í). This subject, introduced 
by Lamb [7], has been pursued with different names (inverse problem, control theory, dinamical 
manipulation, etc. [8-11]). We use dynamical manipulation because somehow it has implicit our 
ideals: to enforce systems to evolve as we want. 

Suppose that a spin state n(í) G S 2 is given. The system of equations (3) for b(í) = /iB(í), 
M(n)b(í) = n(t), is such that detM(n) = 0. Thus, one n(í) does not lead to a unique b(í). If 
b 3 (t) is left arbitrary the other components are: 

6i(í) = [6 3 (í)ni(í) + n 2 (t)]/n 3 (t), b 2 (t) = [b 3 (t)n 2 (t) - n^/n^t). (4) 

As n(í) = i2(í)n(0), b(í) depends of the genèric motion (encoded in R(t) £ SO (3)) and the 
initial condition. This is unsatisfactory because two different fields would be needed to induce 
two trajectòries with common R(t) but different n(0). However, in the direct approach a constant 
field leads unambiguously to R(t). Would it be possible that one R(t) determines a unique b(í) 
when the inverse techniques are applied? This is analysed below for specific forms of R(t). 
Let us suppose first that n(í) rotates with angular velocity 5(t) around 2-axis: 

n(í) = sin 6<o cos(0 o + à(t))i + sin 9 sin(0 o + S(t))j + cos # k, (5) 

where i, j, k are unit vectors along x, y, z, 5(0) = 0. The fields inducing (5) become: 

b(í) = tan 9 [6 3 (í) + 6(t)} [cos(0 o + 6(t))i + sin(0 o + S(t))j] + b 3 (t)k. (6) 

It is clear that b(í) depends on n(0); to avoid this we make b 3 (t) = —S(t). Hence: 

b(í) = -5(t)k. (7) 

If S(t) = ou, we recover the Standard solution for the spin in a constant field. 

Let us analize a more general case including the Standard magnètic resonance model [5] . To this 
aim, suppose that n(í) rotates simultaneously around two fixed directions with variable angular 
velocities. Let us choose one of these directions along k and the other one along a vector e x on 
x — z plane at an angle x from k. The rotation matrix is: 

R(t) = R 3 (P(t))RxHt)) = R3m))R2\-x)R3(-a(t))R 2 (- X ), (8) 

where R 2 (uj) and R 3 (üü) are finite rotation by uj around j and k. If the R 2 (uj) and R 3 (uS) of [1] are 
used to evaluate n(í) = i2(í)n(0) and this is substituted into (4) one arrives at: 

6i(í) = [-v / T^iV2(í)+AiV3]- 1 {sin/3(t){Aà(í)-[63(í)+/3(í)]}iVi(í) 

+ cos^(í) {[6 3 (í) + m][\N 2 (t) + VT^N 3 ] - à(t)N 2 (t)}}, (9) 



b 2 (t) = [-v / T^A2iV 2 (í) + AiV3]" 1 {-cos/?(í){Aà(í)-[Ò3(í)+/j(í)]}iV 1 (í) 

+ sin/3(í) {[6 3 (í) + /3(í)][AiV 2 (í) + v'T^PiVa] - à(t)N 2 (t)}}, (10) 



where À = cosx, n(0) = (xo,yo, Zo), an d: 

iVi(í) = — À sin a(í)xo + cos a(t)y + y/1 — A 2 sin a(t)zo, 

N 2 (t) = A cos a(t)x + sin a(t)y — y/1 — A 2 cos «(t)^, (11) 

iV 3 = VI - A 2 x + Az . 

The field independent of n(0) appears if bs(t) + $(t) = \à(t): 

b(í) = à(í)Vl - A 2 [cos/3(í)i + sin/3(í)j] + [Aà(í) - /3(í)]k. (12) 

There is a consistency condition, n(0) = i?(0)n(0) =>• i2(0) = 1 =>• a(0) = ,3(0) = 0. 

We pass now to the particular cases, firstly those of the introduction. On the one hand, our 
solution (8-12) provides in two ways the case when n(t) rotates around k: by taking x — 0, A = 1, 
and b(í) = [à(t) — $(t)]k; take now a(t) = 0, which leads to b(í) = — $(t)\a. On the other hand, 
if the rotations of n(í) are uniform (a(t) = aot, (3{t) = Pot), we arrive at the traditional model 
used to examine the magnètic resonance [5]: 

b(í) = a Vl - A 2 [cos(/3 í)i + sin(A)í)j] + (Aa - A))k. (13) 

The resonance explanation runs as follows. Suppose that the spin points along k, and it is placed 
in a constant field Bo = B z k = (Acto — A))k/ [i. Hence, at t — the spin state is an eigenstate of the 
'base' Hamiltonian H = — (Aao ~ Po)S z , | -0(0) ) = |+), where S z \+) = (h/2)\+). On S* 2 we have 
n(0) = (0,0, 1). Now, at t = we superimpose to B the field = (a /fx)\/l — \ 2 [cos(f3ot)i + 
sin(/3 í)j]. Thus, a formally 'perturbative Hamiltonian' W = — «oa/I — ^ 2 [cos(f3ot)S x + sm(f3 t)S y } 
is gotten (permitting however the exact treatment). It could induce at t > a transition to the 
orthogonal eigenstate |— ) of H , which on S 2 is represented by the vector n_ = (0,0,-1). The 
probability transition is: 

P+^-(t) = \(-\m)\ 2 = [1 - n 3 (t)}/2 = (1 - A 2 ) [1 - cos(aoí)] /2. (14) 

Notice that >-{t) is small if x i s small or x ~ 7r > i-e-j when the rotations axes of n(í) are 

aligned. The greatest value arises for x = 7r/2 (orthogonal rotation axes): 

P+^_(í) = [l-cos(a í)]/2. (15) 

At r n = (2n + l)7r/ao this probability is 1, i.e., the state certainly will be |— ). This resonance is 
understood analysing the physics in the rotating frame: the rotating observer will see that (due to 
a non-inertial field produced by the rotation) the level spacing between |+) and |— ) will decrease 
with respect to the inertial observer. For A = this spacing will be zero because the non-inertial 
field cancels B , and the transition in the eyes of the rotating observer certainly will be induced. 
In terms of S* 2 there is an equivalent explanation: in the rotating frame it is present a constant 
effective field b e // = ao(y/ï — A 2 i + Ak) around which n(0) = (0,0, 1) starts to precess. But if 
n(0) has to be converted into n_ = (0,0, —1), b e // should have vanishing z component (A = 0). 



In this case, n(í) will simply precess around i with angular velocity ao (period T = 27r/ao), and 
at r n = T/2 + nT = (2n + l)ir/aíQ the spin state will be n(r n ) = n_ = (0, 0, —1). 

Let us remind the absence of essential restrictions on the real functions a(t) and j3(t) in (8-12). 
This makes possible to study a lot of other interesting examples, e.g., the case when the two 
rotations (8) would be induced by a field with a constant third component, 6 3 (í) = bo, which leads 
to a(t) = [b t + P{t))/X. Hence: 



b(í) = ^ [6 + /3(í)][cos/3(í)i + sin/?(í)j] + b k. (16) 

This field has a x — y projection rotating around k with angular velocity $(t) and amplitude 

|&o + /9(í)|^/l/A 2 — 1. It represents a new resonance model analogous to (13-15) where the zero 
Hamiltonian is associated to &ok and the 'perturbation' corresponds to the field component or- 
thogonal to k. Let us suppose that n(0) = (0, 0, 1) and evaluate the probability that at time t the 
system will be in n„ = (0, 0, —1): 

P+-+-(t) = [1 " n 3 (t)]/2 = (1 - A 2 ) [1 - cosa(í)] /2. (17) 

Once again, the greatest probability arises for A = 0: 

P+^_(í) = [1 - cosa(í)]/2, (18) 

and at r n such that a(r n ) = (2n + P H >-(r„) = 1, i.e., the system certainly will make the 

transition from |+) to |— ) due to the 'perturbation' W = — [bi(t)Si + ^(í)^]- 

It is interesting to study the possibility that any n(t) G S 2 would be cyclic, i.e., that the 
evolution operator U{t) (\ip(t)) = U(t)\ip(0))) would be the identity (modulo phase) at some t — r. 
These processes, named evolution loops (EL) [8,11-12], mimic the harmònic oscillator behaviour for 
systems with time-dependent Hamiltonians. The EL are useful to manipulate quantum systems, 
e.g., to induce the squeezing of the wavepacket inside of a Penning trap variant [11]. They are 
helpful as well to produce the rigid displacement of the wavepacket inside a 'magnètic chamber' 
perturbed by homogeneous time-dependent elèctric fields [11]. Some arguments indicate that, 
by applying perturbations, they can be used to induce any unitary operator as the result of the 
precession of the distorted loop [8]. 

For 1/2 spin systems, the EL provide non-trivial cyclic evolutions for which the geomètric phase 
can be explicitly evaluated [13] (see also [12] and below). In our example (8-12), the loop condition 
is immediately formulated: the EL exist if there is a time r for which n(í) simultaneously performs 
n effective rotations around k and / around e x , a(r) = 2ln, (3{t) = 2mr, r > 0, n G Z, / G Z. 
This condition translates into some equations for the field-spin parameters and r. It was proposed 
in [13] for the rotating field (13), and an example of a hypocycloid performing 1 rotation around 
k and 5 rotations around e x was shownF]. Other authors obtained interesting results for the same 
physical system [14]. Here, we want to show 'deformed versions' of the figure in [13], for a field 
of form (16) with f3(t) = aAsin(aot) + (A«i — bo)t. Our results are drawn in figure 1 for two 
essentially different EL with r = 2tt, bo = 3, A = 4/5: a) first we choose a = 0.7, a± = ao = 5; b) 
now we choose a — 1, a± = 5 — 1/27T, ao = 9/4 . 

1 There it was explored also the EL for the spin in the oscillating field b(í) = [bo + òcos(cjí)]k. 
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Fig.l. Spin trajectòries illustrating two EL of period r = 2n induced by (16) with 
(3(t) = a\sm(a t) + (Ac*i - b )t, b = 3, A = 4/5, n(0) = (\/3/2, 0, 1/2), and: a) 
a = 0.7, ai = a = 5; b) a — 1, «1 = 5 — 1/27T, a = 9/4. 

It is easy to see that case a provides a periòdic loop: U(t — 2n) = 1 and U{mr) = 1, m G N. 
This happens because the Hamiltonian, for this choice of parameters, is periòdic with period 
t = 2tï. For case b the parameters are such that U(t — 2ir) = 1, but this EL is not periòdic, 
U(2t) 7^ 1, because the Hamiltonian is not. Different examples of this kind of aperiodic loop, 
when the field oscillates along k, are given elsewhere [13]. 

Now, if an EL is produced at t — r, then any state n(í) becomes cyclic, n(r) = n(0). Thus, 
it makes sense to evaluate the corresponding geomètric phases, which characterize some global 
curvature effects of V = S 2 [3,12-15]. For spin-1/2 systems the geomètric phase turns out to be 
minus one half of the sòlid angle Afl subtended by the oriented closed curve n(í) G S 2 [13], where: 



Aü = 2nn[l - cosxcos(6» - x)} - 2Í7r[l - cos(6» - x)} + sinxsin(6 l - x) pcosadt, (20) 



where we have taken n(0) = (sin#, 0, cos 8). This expression is equal to the corresponding one for 
the rotating field (13) but for the last term [13], which depends on the explicit form of both a(t) 
and f3{t). If 63 (í) = 60 (the case discussed at (16-18)) we get: 



AÍ7 = 2r7,7r[l — cosxcos(6 ) — %)] — 2/7r[l — cos(6 ) — x)] — 6 sinxsin(6' — x) / cosadí, (21) 




(19) 



For our case (8-12), this formula can be evaluated using the loop condition: 



■T 




and the integral depends just of a(t). Notice that this integral vanishes if a(t) = a t and a r = 
21it, which leads once again to the known result for the EL in the Standard case (13). The simplest 



situation leading to EL for which the integral in (21) does not vanish arises when a(t) is quadratic 
in t. With this choice, that term will involve Fresnel integrals^. 

In conclusion, the manipulation techniques are appropriate to generate solvable examples for 
the spin evolution in homogeneous magnètic fields. When these techniques are applied to produce 
the EL, it arises cyclic evolutions for which the geomètric phases can be explicitly evaluated. 
We hope that our treatment has shed as well some light on the functioning of the resonance 
mechanism. 
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